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Averaging Functions 
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N. Santiago^" 



Abstract 

In this paper we propose a special type of aggregation function which generalizes the notion of Ordered 
Weighted Averaging Function - OWA. The resulting functions are called Dynamic Ordered Weighted Averaging 
Functions — DYOWAs. This generalization will be developed in such way that the weight vectors are variables 
depending on the input vector. Particularly, this operators generalize the aggregation functions: Minimum, Maximum, 
Arithmetic Mean, Median etc, which are extensively used in image processing. In this field of research two problems 
are considered: The determination of methods to reduce images and the construction of techniques which provide 
noise reduction. The operators described here are able to be used in both cases. In terms of image reduction we 
apply the methodology provided in [1], We use the noise reduction operators obtained here to treat the images 
obtained in the first part of the paper, thus obtaining images with better quality. 

Index Terms 

Aggregation functions, OWA functions. Image reduction. Noise reduction 



I. Introduction 

Image processing has great applicability in several areas. In medicine, for example, they can be applied 
to: Identify tumors [2]; support techniques in advancing dental treatments [3], etc. Such images are not 
always obtained with a suitable quality, and to detect the desired information, various methods have been 
developed in order to eliminate most of the noise contained in these images. 

Another problem addressed in image processing is the decrease of resolution, usually aiming the 
reduction of memory consumption required for its storage [4], 

* + 
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There are several techniques for image reduction in the literature, more recently Paternain et. at. [1] 
constructed reduction operators using weighted averaging aggregation functions. The proposed method 
consists of: (1) To reduce a given image by using a reduction operator; (2) To build a new image from 
the reduced one, and (3) To analyze the quality of the last image by using the measures PSNR and MSIM 
[4]. 

In this work we introduce a class of aggregation functions called: Dynamic Ordered Weighted 
Averaging Function - (DYOWA). They generalize the OWA function introduced by Yager [5] and in 
particular the operators: Arithmetic Mean, Median, Maximum, Minimum and cOWA. We provide a range 
of their properties such as: idempotence, symmetry and homogeneity as well two methods 1 : (1) for image 
reduction and (2) for noise treatment. 

This paper is structured in the following way: Section 2 provides some basics of Aggregation Functions 
Theory. SECTION 3 introduces Dynamic Ordered Weighted Averaging functions, shows some examples 
and properties, and introduces a particular DYOWA function, called H. which will be fundamental for 
this work. In Section 4 we provide an application of DYOWA ' s to image reduction and in Section 
5, we show that DYOWA functions are able to treat images with noise, aiming to improve the reduction 
method used in Section 4. Finally, section Section 6 gives the final remarks of this work. 

II. Aggregation Functions 

Aggregation functions are important mathematical tools for applications in several fields: Information 
fuzzy [6]; Decision making [7]-[ll]; Image processing [1], [2], [12] and Engineering [13]. In this section 
we introduce them together with examples and properties. We also present a special family of aggregation 
functions called Ordered Weighted Averaging - OWA and show some of its features. 

A. Definition and Examples 

Aggregation functions associate each entry x with n arguments in the closed interval [0, 1] an output 
value also in the interval [0, 1]; formally we have: 

Definition 1. An n-ary aggregation function is a mapping f : [0, l] n — > [0, 1], which associates each 
n-dimensional vector x to a single value /(x) in the inten’al [0, 1] such that: 

1) /(0, ..., 0) = 0 and /(l, ..., 1) = 1; 

‘These methods were implemented by using Java 1.8.0_31 software in a 64 bits MS Windows machine. 
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2) If x < y, i.e., Xi < y it for all i — 1, 2, n, then /(x) < /( y). 

Example 1. 

(a) Arithmetic Mean: Arithix) = — (x\ + x 2 ... + x n ) 

n 

(b) Minimum: Mm(x) — min{x i, X 2 , 

(c) Maximum: Max(x ) = max{x\ } x 2 , ..., x n }; 

n 



(d) Harmonic mean: / n (x) = — y- 

X\ X2 Xn 

From now on we will use the short term “aggregation” instead of “n-ary aggregation function”. 
Aggregations can be divided into four distinct classes: Averaging, Conjunctive, Disjunctive and Mixed. 
Since this paper focus on averaging aggregations, we will define only this class. A wider approach in 
aggregation can be found in [14]— [17]. 

Definition 2. An aggregation is called Averaging, if for all x 6 [0, l] n we have: 



Mm(x) < /(x) < Max(x) 



Example 2. The Arithmetic Mean, the Maximum and the Minimum are averaging aggregation functions. 



B. Special Types Aggregation Functions 
An aggregation function /: 

(1) is Idempotent if, and only if, f(x, ...,x) = x for all x G [0, 1]. 

(2) is Homogeneous of order k if, and only if, for all A e [0, 1] and x e [0, l] n , f(Xxi , Ax 2 , ..., Xx n ) = 

X k f(xi,x 2 , ... ,x n ). When / is homogeneous of order 1 we simply say that / is homogeneous. 

(3) is Shift-invariant if, and only if, f(x\ + r,x 2 + r, .. ,x n + r ) = f(x i,x 2 , ..,x n )+r, for all r e [—1, 1], 
x G [0, l] n such that (xi +r,x 2 + r, ..., x n + r) G [0, l] n and f(x i,x 2 , ..., x n ) + r G [0, 1]. 

(4) is Monotonic if, and only if, x < y implies /(x) < /( y). 

(5) is Strictly Monotone if, and only if, /(x) < /( y) whenever x < y, i.e. x < y and x f y. 

(6) has a Neutral Element e G [0, 1], if for all t G [0, 1] at any coordinate input vector x, it has to be: 

/(e, ..., e, t, e, ... ,e ) = t, and 



(7) / is Symmetric if, and only if, its value is not changed under the permutations of the coordinates 
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of x, i.e, we have: 



f(x i,x 2 , ...,x n ) = f 0 




For all x and any permutation P : {1, 2, n} — > {1, 2, n}. 

(8) An Absorbing Element (Annihilator) of an aggregation function /, is an element a G [0, 1] such 
that: 



(9) A Zero Divisor of an aggregation function is an element a 6 ]0, 1[, such that, there is some vector 
x with Xj > 0, for all 1 < j < n, and f(x i, Xj_ i, a, x 3+l , x n ) = 0. 

(10) A One Divisor of an aggregation function / is an element a € ]0, 1[ such that, there is some vector 
x with Xj < 1 , for all 1 < j < n, and f(x i, Xj-i, a, x 3+ \, x n ) = 1 . 

(11) If N : [0,1] — >■ [0,1] is a strong negation 2 and / : [0, l] n — >■ [0,1] is an aggregation function, then 
the dual aggregation function of / is: 



which is also an aggregation function. 

Example 3. 

(i) The functions: Arith, Min and Max are examples of idempotent, homogeneous, shift-invariant, 
monotonic and symmetric functions. 



ladors. 

(iii) Min, Max and Arith do not have zero divisors and one divisors. 

(iv) The dual of Max with respect to negation N(x) = 1 — x is the Min function. 

C. Ordered Weighted Averaging Function - OWA 

In the held of aggregation functions there is a very important subclass in which the elements are 
parametric; they are called: Ordered Weighted Averaging or simply OWA [5]. 

An OWA is an aggregation function which associates weights to all components x t of an input vector 



f (:i i, X 2 , ■•*, Xi—i, a , Xi~\~i, •*., Xr,fj a 



f d (x 1 ,x 2 , x n ) = N(f(N(x 1 ),N(x 2 ), ..., N(x n ))) 



(ii) Min and Max have 0 and 1 elements as annihilator, respectively, but Arith does not have annihi- 



x. To achieve that observe the following definition. 

2 A strong negation is an antitonic function N : [0, 1] — » [0, 1] such that N(N(a)) = a for all a G [0, 1]. 
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Definition 3. Let be an input vector x = (x\,X 2 ,...,x n ) G [0, l] n and a vector of weights w = 

n 

(u>i, . . . ,w n ), such that = 1. Assuming the permutation: 

i = 1 

Sort(pt) ( 2-^(1) ; 3'p(2) j ■ ■ ■ )2Jp(n)) 



such that x p (i) > x p ( i+1 ), i.e., x p ^ > x p ( 2 ) > ••• > x p ( n y the Ordered Weighted Averaging (OWA) 
Function with respect to w, is the function OWA w : [0, l] n -4 [0, 1] such that: 



OVFA w (x) = £>■ x p {i) 

i=l 

In what follows we remove w from OWA w (x). The main properties of such functions are: 

(a) For any vector of weights w, the function OWA(x) is idempotent and monotonic. Moreover, OWA(x) 
is strictly increasing if all weights w are positive; 

(b) The dual of a OWA w is denoted by (OWA) d , with the vector of weights dually ordered, i.e. 
(■ OWA w ) d = OWA w d, where w d = (w p{n) ,w p{n _ ...,w p{1) ). 

(c) OWA are continuous, symmetric and shift-invariant functions; 

(d) They do not have neutral or absorption elements, except in the special case of functions OWA of 

Max and Min. 

1 ) Examples of special functions OWA: 

1. If all weight vector components are equal to -, then OWA(x) = Arith((x). 

2. If w = (1, 0, 0, ..., 0), then OWA(x) = Max(x ). 

3. If w = (0, 0, 0, ..., 1), then OWA(x) = Mm(x). 

4. if Wi = 0, for all i, with the exception of a k - th member, i.e, u;/,. = 1, then this OWA is called static 

and OWA w (x) = Xk 

5. Given a vector x and its ordered permutation Sort(x) = (x(i), ■ ■ ■ ,X(n)), the Median function 



Med(x) 



|(x (fc) + x (fc+ i)), if n = 2k 
X( k + 1), ifn = 2fc + l 



is an OWA function in which the vector of weights is defined by: 

• If n is odd, then Wi = 0 for all i f |~|~| and w \ n / 2 ] = 1. 

• If n is even, then = 0 for all i f and i f L^J’ an d w [n/2] = w\n/2\ = \- 



Example 4. The n-dimensional cOW A function [18] is the OWA operator, with weighted vector defined 
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by: 

• If n is even, then Wj = for 1 < j < and w n / 2+i = w n / 2 -i+i, for 1 < i < |. 

. If n is odd, then Wj = 2(2 ) £ i> . for 1 < j < w n / 2+ i = w n / 2 -i+i, for\<i<\, and W( n +i )/2 = 

(n— 1)/2 

1-2 £ Wi . 

j = i 

The OWA functions are defined in terms of a predetermined vector of weights. In the next section we 
propose the generalization of the concept of OWA in order to relax the vector of weights. To achieve that 
we replace the vector of weights by a family of functions. The resulting functions are called Dynamic 
Ordered Weighted Avegaring Functions or in short: DYOWAs. 

III. Dynamic Ordered Weighted Avegaring Functions - DYOWA 
Before defining the notion of DY OWA functions, we need to establish the notion of weight-function. 
Definition 4. A finite family of functions Y = {/,; : [0, l] n — » [0, 1] | 1 < i < n} such that: 

n 

5^/i( x ) = i- 

i=l 

is called family of weight-function (FWF). 

A Dynamic Ordered Weighted Averaging Function or simply DYOWA associated to a FWF T is a 
function of the form: 

n 

OyWA r (x) = ^/ i (x)-T, 

i= 1 

Below we show some examples of DYOWA operators with their respective weight-functions. 

Example 5. Let beT = {/i(x) = 1 | 1 < i < n}. The DY r OWA operator associated to T, DYOWAy^x), 
is Arith{x). 

Example 6. The function Minimum can be obtained from V = { f , \ 1 < i < n}, where /i(x) = f 2 {x ) = 

■ • • = fn—i (x) = 0 and / n (x) = 1, for all x G [0, l] n . 

Example 7. Similarly, the function Maximum is also of type DYOWA with Y dually defined. 

Example 8. For any vector of weights w = (wi, w 2 , ..., w n ), A function 0I17l w fx) is a DYOWA in 
which the weight -functions are given by: /, ; (x) = w p ^, where p : {1,2 , — > {1, 2, ■ ■ ■ , n} is 

the permutation, such that p{i) = j with x, = x^y For example: If w = (0.3, 0.4, 0.3), then for x = 
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(0.1, 1.0, 0.9) we have x\ = x ( 3 ), £2 = £(i) and £3 = £( 2 ). Thus, /i(x) = 0.3, / 2 (x) = 0.3, /s(x) = 0.4, 
and DYOWA(Y) = 0.3 • 0.1 + 0.3 • 1.0 + 0.4 • 0.9 = 0.69 

Remark 1. Example 8 shows that the functions OWA, introduced by Yager, are special cases of DYOWA 
functions. There are, however, some DYOWA functions which are not OWA. 

Example 9. Let V = {sin(£) ■ y, 1 — sin(£) • y}. The respective DYOWA function is DYOWA{x, y ) = 
(sin(£) ■ y) ■ x + (1 — sin(£) • y) ■ y, which is not an OWA function. 



A. Properties of DY OWA Functions 

The next theorem characterizes the DYOWA functions which are also aggregations. 

Theorem 1. Let V = {/ 1 , • • • , f n } be a FWF. A DYOWAy is an aggregation function if, and only if, it 
is monotonic. 

Proof. Obviously, if DY'OWAy is an aggregation, then it is monotonic function. Conversely, if DYOWAy 
is monotonic, then for it to become an aggregation, enough to show that 

DYOWA r (0, ..., 0) = 0 e DYOWA r ( 1, 1) = 1, 

this follows from the definition of DY'OWA. □ 

Corollary 1. A DY OWA is an aggregation function if and only if, it is an a aggregation of type averaging. 
Proof. For all x = (x-y, ... ,£ n ) have to 

Mm(x) < Xi < Max(x), \/i = 1,2 ,...,n. 



So, 



^2 /*( x ) ■ /<( x ) • x i < ^2 -ft( x ) ■ 

2=1 2=1 2=1 



n 

but as /*( x ) — 1. it follows that 

2=1 



n 

Mm(x) < ^ /j(x) • Xj < Max(x ) 
2=1 



□ 
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Corollary 2. All functions of the type DY OWA presented in examples 4 , 5, 6, 7 and 8 are averaging 
aggregation functions. 

Proof Just see that those functions are monotonic. □ 

Proposition 1. For every T, DYOWAy is idempotent. 

Proof If x = (x, ..., x) with t G [0, 1], then: 

n n 

DYOWArfx) = ^ /-(x) • x = x • /j(x) = x 

2=1 2=1 

□ 



This property is important because it tells us that every DYOWA is idempotent, regardless it is an 
aggregation or not. 

Proposition 2. IfT is invariant under translations, i.e, f(x i + A, X 2 + A, ...,x n + A) = ffx i,X 2 , x n ) 
for any x G [0, l] n , for i G {1.2,---, n} and A G [—1, 1], then DYOWAy is shift-invariant. 

Proof. Let x = (xi, x n ) G [0, l] n and A G [—1, 1] such that (xy + A, x 2 + A, x n + A) G [0, l] n . then, 

DYOWAy {x,\ + A, x n + A) = 

n 

— fi {% i + A, x n + A) • ( Xi + A) 

2=1 

n 

= Y. fi( x i + A, X n + A) • Xi 

2=1 

n 

+ fi(% i + A, x n + A) • A 

2=1 

n 

= -Xj + A 

2=1 

DYOWA r (xy, ..., x n ) + A 



□ 



Proposition 3. If Y is homogeneous of order k (i.e., if f is homogeneus of order k, for each f E T), 
then DY OWAy(x) is homogeneous of order k + 1. 
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Proof. Of course that, if A = 0, then DYOWA r (Xxi , ..., \x n ) = Xf(x\, ...,x n ). Now, to A 7^ 0 we have: 

n 

DYOWA r (\xi,...,\x n ) = y^ffXxi, ...,Xx n ) ■ Xxj 



i= 1 



A • ^ X k f i (x 1 , ...,x n )xi 



i=l 



= X k+1 ■ DYOWA r ( Xl , ..., x n ) 



□ 



Example 10. Let V be defined by 
Then, 



E x 3 
3 = 1 



if x 1 = ... = x n = 0 
- , otherwise 

if x 1 , . x n = 0 



DYOWAr(x) = 



0 , 



E x i 

Ar — , otherwise 
E x i 

i = 1 



This DYOWAy is idempotent, homogeneous and shift -invariant. However, DY OW Ay is not monotonic, 
since DYOWA r (0.5, 0.2, 0.1) = 0.375 and DYOWA r (0.5, 0.22, 0.2) = 0.368. 

The next definition provides a special FWF, which will be used to build a DYOWA whose properties 
are very important for this paper. 



Definition 5. Consider the family I of functions 



If x = (x, ...,x) 



/*(x) = < 



— —r [ 1 TT- 1 I , otherwnse 

n 1 E |xj— Afed(x)| 1 

3 = 1 



Then, T is a FWF, i.e. / ( (x) = 1, for all x e [0, l] n . Let H /?c //n? associated DYOWA. The computation 
2=1 

of FI can he performed using the following expressions: 



if x = (a, ...,:r) 



H(x) = 



a;* - Med ( x )l , otherwise 

1 E |*i— Med(x)| 1 

J = 1 
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Example 11. Let be n = 3. So, for x = (0.1, 0.3, 0) we have 

/i(x) = 0.5, ^ 2 (x) = 0.167, fs(x) = 0.333 and H(x) = 0.08. 



Proposition 4. The weight-functions defined in Definition 5 are such that: j,(x\ + A, .... x n + A) = 
fi(x i,x 2 ,...,x n ) and ffiXxi, Xx n ) = ffixx , ..., x n ), for any 1 < i < n. 



Proof Writing x' = (xi+A, ..., x n +X), then f(x±+X , ..., x n +X) = (/i(x'), ..., / n (x')). Clearly, Medfx!) = 
Med{x) + A. Thus, for x f (x, ..., x) we have: 



m *) = ^ h- 



|®i+A— Med{x!) | 

n 

53 \xj+\—Med(-x!)\ 
i = l 



1 

n— 1 




\xj+\—(Med(x)+ A) | 

n 

Y \xj-\-X—(Med(x.)-\-X)\ 
3 = i 



1 

72—1 




\xi—Med{x)\ 

n 

Y | Xj—Med{x) 

3 = 1 



= /*(x). 



Therefore, /(x') = (fi(x'), ..., / n (x')) = (/i(x), ..., /„(x)). The case in which x = (x, ..., x) is immediate. 

To check the second property, make x" = (Axi, ..., Ax n ), note that Med(x") = Xined(x) and for 
X ^ (x, ...,x) 



MY = Mr 



1 | , \\xi~ Med(Ax)| 

53 Axj —M ed(Ax) 
i=l 



1 | i |A:Ci— AMed(x)| 

n — i I ^ 7™ 

53 |Axj— AMed(x)| 

3=1 



1 | 1 _ |A[-|xj— Med(x)| 

r7 — 1 I 1 n 

|A|- 53 \xj-Med(x) 
3 = 1 



1 I i _ |xj-Med(x)| 
n_1 53 |xj— Med(x)| 

= /*« 



Therefore, /(x") = (/i(x"), ..., /„(x")) = (/i(x), ..., / n (x)) = /(x). The case in which x = (x,...,x) is 
also immediately □ 
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li 



Corollary 3. H is shift-invariant and homogeneous. 



□ 



Proof. Straightforward for propositions 2 and 3. 

The function H is of great importance to this work, since this function, as well as some DYOWA’s 
already mentioned will provide us tools able: (1) To reduce the size of images and (2) To deal with noise 
reduction. 

Now, we present some other properties of function H. 

B. Properties of H 

In addition to idempotency, homogeneity and shift-invariance H has the following proprerties. 
Proposition 5. H has no neutral element. 

Proof. Suppose H has a neutral element e, find the vector of weight for x = (e, x,e , ... , e ). Note 

that if n > 3, then Med(x) = e and therefore, 



x = 



1 | -i _ \xj—Med(x)\ 

77 — 1 I n 

E |®j— Afed(x) | 

3 = 1 



n—1 



l - 



\xj-e\ 

n 

E \xj-e | 

3 = 1 



therefore, 



i.e.. 



and 



_ _!_ (\ _ h»~ e l 

n—1 \ \x— e\ 



, — -r, if Xi — e 

/i(x) = { , to n > 3 

0, if Xi = x 



fix') = (A_ n 

AN Vn-l > ■■■> n-1 ’ n-1 ’ n-l) 



H(x) = (n — 1) 



n — 1 



= e 



But since e is a neutral element of H, H(x) = x. Absurd, since we can always take i^e. 

x c 

For n = 2, we have Med(x) = , where x = (x, e) or x = (e, x). In both cases it is not difficult 

2 

x c 

to show that /(x) = (0.5, 0.5) and H(x) = — - — . Thus, taking x f e, again we have H(x, e) x. □ 
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Proposition 6 . H has no absorbing elements. 

X\ 4“ X2 

Proof. To n = 2, we have H(x) = — - — , which has no absorbing elements. Now for n > 3 we have 
to x = (a, 0 , ..., 0 ) with Med(x) = 0 therefore, 

1 



1 - -) = 0 and fi = — — -,Vi = 2 ,...,n. 
n — IV a / n — 1 



A( x ) = T (i-“) 

n — IV aJ 



therefore, 



H(a, 0,...,0) = 0 • a 4 • 0 + ... 4 • 0 = a =4* a = 0, 



n — 1 



n — 1 



but to x = (a, 1, ..., 1) we have to Med(x) = 1. Furthermore, 

fi (x) = — — (l - . - a ) = 0 



n — 1 



and 



fi = 



n — 1 



para z = 2,3, ..., n. 



therefore, 



H(a, 1 , . .. , l) — 0 • 0 4- , • 1 T • • • 4~ ~ ■ 1 — a a — 1 . 



n — 1 



n — 1 



With this we prove that H does note have annihiladors. 



□ 



Proposition 7. H has no zero divisors. 

n 

Proof. Let a E ]0, 1[ and consider x = (a, x 2 , x n ) e ]0, l] n . In order to have H(x) = ffx.) ■ Xi = 0 

i= 1 

we have /i(x) • Xi = 0 for all i = 1,2, ...,n. But as a 7 ^ 0 and we can always take x 2 , x 3 , ...,x n also 
different from zero, then for each i — 1,2, .... n there remains only the possibility of terms: 



fi (x) = 0 para i = 1 , 2 , ...,n. 

n 

This is absurd, for /i(x) e [0, 1] e /i(x) = 1. like this, H has no zero divisors. □ 

2=1 

Proposition 8. H does not have one divisors 

Proof. Just to see that a e ]0, 1[, we have to H(a, 0, ..., 0) = /i(x).a < a < 1. □ 

Proposition 9. H is symmetric. 

Proof. Let P : {1, 2, n} — > {1, 2, ..., n} be a permutation. So we can easily see that Med(x P i y i),x P ( 2 ), ..., xp( n )) = 
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Med(x i, x 2 , ...,x n ) for all x = (xi, x 2 , ■■■, x n ) E [0, l] n . We also have to \x P ^-Med(x P ^), x P(2 ), ...,x P ( n ))| 



i = 1 



^2 | Xi — Med(x) |. Thus, it suffices to consider the case where (x P (i), x P ( 2 ), ..., x P ( n )) f (x,x, ...,x). But 



2—1 



(x P (i), x P(2 ), x P(n) ) f (x, x , x) we have to: 

H(x P (i), X P ( 2 ), ..., X P ( n )) = 

1 I ™ Xp (i) \xp (i) -Med(xp (1) ,...,x p(n) )\ 

~ n-1 2^ X P(i) 



2—1 



E X P(i) 

i = 1 



n — 1 


M 
1 e 

1 


n 

E x 4 

i= 1 


1 'v 


n — 1 


n-1 ^ 

i- 


n 

E x i 

i= 1 


i \ 



E \ x P(i)~ Med ( x P(i),-i x P(n))\ 
3 = 1 



E x P ^\xp^-Med(xi,...,x n )\ 
n 

2—1 E |^p(i)-Med(^i,...,tr n )| 

j = i 



n 

E \xi-Med(xi,...,x n )\ 

3 = 1 

a; i |xi—ALed(a:i, ...,a; n ) | 

n 

2 — 1 E |®»— Med(®i,...,;E n )| 
j=i 



= H(xi, ...^ 



□ 



Therefore, H satisfies the following properties: 

. Idempotence 
. Homogeneity 

• Shift-invariance 

• Symmetry. 

• H has no neutral element 

• H has no absorbing elements 

• H has no zero divisors 

. H does not have one divisors 

Remark 2 . Unfortunately we do not prove here the monotonicity of H, due to its complexity, but we 
suspect that it is true. This demonstration will be relegated to a future work. 

The next two sections show the suitability of DY OWA. They will provide applications for image and 



noise reduction. 
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IV. DYOWA’S AS IMAGES REDUCTION TOOLS 



In this part of our work we use the functions DYOWA studied in Examples 4, 5, 6, 7 and 8, and 
definition 5 to build image reduction operators, the resulting images will be compared with the reduced 
image obtained from the operator function H. 

An image is a matrix mxn, M = A (i , j ), where each A(i, j ) represents a pixel. In essence, a reduction 
operator reduces a given image m x n to another m! x n', such that m! < m and n' < n. For example, 

0.1 0.2 0 0.5 

0.3 0.3 0.2 0.8 
1 0.5 0.6 0.4 

0 0.3 0.5 0.7 

In Grayscale images the value of pixels belong to the set [0, 255], which can be normalized by dividing 
them by 255, so that we can think of pixels as values in the range [0, 1]. 



0.1 0 
1 0.6 




Fig. 1. Example of image in Grayscale. 



There are several possible ways to reduce a given image, as shown in the following example: 

Example 12. The image 



0.8 


0.7 


0.2 


1 


0.5 


0.5 


0.6 


0.2 


0.3 


0.1 


1 


0 


0 


0 


0.6 


0.4 


0.9 


1 


0.1 


0.2 


0.3 


0.4 


0.5 


0.6 



Page 15 of 24 IEEE Transactions on Fuzzy Systems Proof 

JOURNAL OF RTeX CLASS FILES, VOL. , NO. , 2015 



15 



1 

2 

3 



14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 

59 

60 



can be reduced to another 2x3 by partitioning M in blocks 2 x 2: 



4 

5 
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M = 














10 

11 






0 


0 
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12 
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and applying to each block, for example, the function fix, y. z, w) = Max(x,y, z,w): We obtain, the 
image: 

M* = 



0.8 1 1 
0.2 0.6 1 



Applying g(x,y, z,w) = Min(x,y, z,w) we would obtain: 



M** = 



0.2 0.1 0 

0 0.2 0.5 



In fact, if we apply any other function, we get a new image (usually different from the previous one), 
but what is the best? 

One possible answer to this question involves a method called magnification or extension (see [19]— 
[21]), which is a method which magnifies the reduced image to another with the same size of the original 
one. The magnified image is then compared with the original input image. 

Example 13. From M* and M**, we get images 4x6, M' and M", simply cloning each pixel , 

x x 
x x 

We obtain, new images 



M 1 = 
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